Final results from the Palo Verde Neutrino Oscillation Experiment 
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The analysis and results are presented from the complete data set recorded at Palo Verde between 
September 1998 and July 2000. In the experiment, the Uc interaction rate has been measured at 
a distance of 750 and 890 m from the reactors of the Palo Verde Nuclear Generating Station for a 
total of 350 days, including 108 days with one of the three reactors off for refueling. Backgrounds 
were determined by (a) the swap technique based on the difference between signal and background 
under reversal of the positron and neutron parts of the correlated event and (b) making use of the 
conventional reactor-on and reactor-off cycles. There is no evidence for neutrino oscillation and the 
mode z/c — > Vx was excluded at 90% CL for Am^ > 1.1 x 10""^ eV^ at full mixing, and sin^ 29 > 0.17 
at large Am^. 



I. INTRODUCTION 

We report here the final results of a long baseline study 
of i>c oscillations at the Palo Verde Nuclear Generating 
Station. This is a continuation of the work reported ear- 
lier in ||lj,y in which details of the experiment and first re- 
sults are described. Hence we only briefly describe the de- 
tector and the analysis, stressing the improvements and 
final results. Since the previous report, the data sample 
has been almost doubled. Improvements have been made 
on reconstruction and simulation, reducing the system- 
atic error by one-third. 

The experiment was originally motivated by the obser- 
vation of an anomalous atmospheric neutrino ratio v^/vc 
reported in several independent experiments P-p[. The 
mass parameter suggested by this anomaly is in the range 
of lO^'^ < Ato^ < 10^'^ eV^ for two-flavor neutrino os- 
cillation. The Palo Verde experiment, together with the 
Chooz experiment 0,^] with a similar baseline, were able 
to exclude v^ — > v^ oscillations as the dominant mecha- 
nism for the atmospheric neutrino anomaly. While the 
experiment has pursued its goal of exploring the then 
unknown region of small Attt,^ , recent data from Super- 
Kamiokande iHl favor the v., -^ v^ oscillation channel 



II. EXPERIMENTAL TECHNIQUE 

A. Detector and signal 

The experiment was performed at the Palo Verde Nu- 
clear Generating Station in Arizona. The plant consists 
of three identical pressurized water reactors with a total 
thermal power of 11.63 GW. The detector was located 



at a shallow underground site, 890 m from two of the 
reactors and 750 m from the third. The 32 meter-water- 
equivalent overburden entirely eliminated any hadronic 
component of cosmic radiation and reduced the cosmic 
muon flux. 

The segmented detector Q consisted of 66 acrylic cells 
filled with 11.34 tons of Gd-loaded Hquid scintillator [|j. 
A 0.8 m long oil buffer at both ends of each 9 m long cell 
and aim buffer filled with water (105 t) surrounding the 
central detector shielded it from radioactivity originat- 
ing in the photomultiplier tubes (PMTs) and laboratory 
walls as well as from neutrons produced by cosmic muons 
passing outside of the detector. The outermost layer of 
the detector was an active muon veto counter, providing 
47r coverage. All materials used in the construction of the 
detector and the laboratory were screened for there ra- 
dioactivity content by means of low background gamma 
ray spectroscopy in order to control detector background. 

The t'e flux was detected via the correlated positron 
and neutron sub-events from the reaction PeP^ne"'". The 
sub-events are (1) the positron's kinetic energy (< E >~ 
2.4 MeV) and two prompt annihilation 7's and (2) the 
subsequent (delayed with a time constant of ^ 27/is ) 
7's from capture of the thcrmalized neutron on Gd (with 
energy ^ 8 MeV). 

The data acquisition (DAQ) electronics was built as 
a dual bank system, allowing both parts of the sequen- 
tial Vc capture event to be recorded with no deadtime 
by switching between banks. Signals from each PMT 
were discriminated by two thresholds: a high threshold 
corresponding to '^600 keV for energy deposits in the 
middle of the cell and a low threshold corresponding to 
^--^40 keV also in the middle of the cell, or a single photo- 
electron (SPE) at the PMT. The trigger processor, a field 
programmable gate array, searched for triple patterns in 



the central detector for each of the sub-events, requiring 
one high and at least two low discriminator signals from 
neighboring cells |1^. All events with two triple signals 
within 450 (is of each other were written to disk. 

A veto signal following the passage of a muon (typi- 
cal veto rates were ~2 kHz) disabled the central detec- 
tor trigger for 10 /iS. With each event, the time and 
hit pattern of the previous muon in the veto counter 
was recorded for offline use along with information as 
to whether or not the muon passed through the target 
cells. The veto inefficiency was measured to be 2.5±0.2% 
for stopping muons (one hit missed) and 0.07±0.02% for 
through-going muons (two hits missed). 

In order to maintain constant data quality during run- 
ning, a protocol of continuous calibration and monitoring 
of all central detector cells was followed. Relative tim- 
ing and position were calibrated with blue light emitting 
diodes (LEDs) installed inside each cell. Additional blue 
LEDs illuminated optical fibers at the end of each cell, 
providing information about PMT linearity and short 
term gain changes. LED and fiber optic scans were per- 
formed once a week. For absolute energy calibration and 
determination of the positron and neutron detection ef- 
ficiencies, as well as mapping of the light attenuation 
in each cell, radioactive sources were used. A complete 
source scan was undertaken every 2-3 months. Further 
details have been described in Ref. ||]. 

B. Expected i'c interaction rate 

To evaluate the expected Dc interaction rate in the de- 
tector, the power and fuel composition of the three re- 
actors must be known. The calorimetric methods, based 
on the measurement of temperature and water flow-rate 
in the secondary cooling loop, provided power determi- 
nation with 0.7% uncertainty. 

The fission rates in the three reactor cores were calcu- 
lated daily using a simulation code provided by the man- 
ufacturer of the reactors. The output of the core simula- 
tion had been checked by measuring isotopic abundances 
in expended fuel elements in the core; errors in fuel ex- 
posure and isotopic abundances are estimated to cause 
< 0.3% uncertainty in the Po flux estimate. Four isotopes 
— 239pu, 24ipu, 235u, and 238U — producc virtually aU 
the thermal power as well as all the i>o's. Measurements 
of the neutrino yield per fission and energy spectra exist 
for the first three isotopes 0,|2|. The '^^^V yield, which 
contributes 11% to the final j/o rate, was calculated from 
theory p3| with an uncertainty of 10%. The contribu- 
tion of 238^ -j-Q ^j^g uncertainty of the total neutrino rate 
is therefore ~1%. 

The Do energy spectrum was reconstructed from the 
measured positron kinetic energy. The approximate re- 
lation Ejy^ = E^+ -f 1.8 MeV is slightly modified by the 
kinetic energy carried away by the neutron (~50 keV). 



The cross section of the detection reaction is accurately 
known p4| ; the dominant uncertainty (0.2%) stems from 
the neutron lifetime. 

Previous short baseline reactor experiments have found 
good agreement between calculated and observed neu- 
trino fluxes p5|-p^. In particular, Ref. [0 quoted an 
uncertainty in the neutrino flux per fission of 1.4%. To- 
gether with the combined uncertainty of 1.5% of the re- 
actor power, the distance to the detector and the number 
of target atoms, the total systematic uncertainty of the 
De interaction rate therefore amounts to 2.1%. 

The expected v^ interaction rate in the whole target, 
both scintillator and the acrylic cells, is plotted in Fig. |l| 
for the case of no oscillation from July 1998 to July 2000. 
Around 220 interactions per day are expected with all 
three units at full power. Four periods of sharply re- 
duced rate occurred when one of the three reactors was 
off for refueling, the more distant reactors contributing 
each approximately 30% of the rate and the closer reactor 
the remaining 40%. The short spikes of decreased rate 
are due to accidental reactor outages, usually less than 
a day. The gradual decline in rate between refuelings 
is caused by fuel burn-up, which changes the fuel com- 
position in the core and the relative fission rates of the 
isotopes, thereby affecting slightly the yield and spectral 
shape of the emitted Dc fiux. 
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FIG. 1. The calculated i>c interaction rate in the detector 
target for the case of no oscillations. The four long periods of 
reduced flux from reactor refuelings were used for background 
subtraction. The decreasing rate during the full power opera- 
tion is a result of the changing core composition as the reactor 
fuel is burned. 



III. MONTE CARLO SIMULATION 
A. Detection efficiency 

Accurate understanding of the v^ efBciency is crucial. 
Therefore, as described in y], two parallel and indepen- 
dent event reconstruction and detector simulation Monte 
Carlo codes have been developed. Both give consistent 
results; in |g] the reported results were based on one of 
them, here most of the results are based on the second 
method. 

A Monte Carlo model with a detailed simulation of 
the detector response, including the PMT pulse shape, is 
essential to simulate the rather strong dependence of the 
i>c efficiency on the event location in the detector and, 
to a lesser extent, on time due to some scintillator aging. 
A variety of measurements was performed to cross check 
the Monte Carlo modeling of the detector response. 

The simulation code pi contains the whole detector 
geometry and simulates the energy, time, and position of 
energy deposits in the detector using GEANT 3.21 [ pi[ . 
GFLUKA [E2| is used to simulate hadronic interactions, 
while for the low energy neutron transport GCALOR |2^] 
is employed. Scintillator light quenching, parameterized 
as a function of ionization density, is included in the sim- 
ulation J24|. 

Given the output of the physics generators, the Monte 
Carlo simulates the detector response in the form of PMT 
pulses which arc converted into time and amplitude digi- 
tizations and trigger hits. Digitized data are then recon- 
structed with the same programs as real data, providing 
the trigger and selection cuts efficiencies. 



B. Improvements of the simulation 

Since the initial results were published g, the data 
sample has been almost doubled. There were also im- 
provements in the analysis due to refinements of the sim- 
ulation of the detector response. Three changes had the 
largest impact on the quality of the simulation: 

• The pulse-shapes of several hundred single- 
photoelectron (SPE) signals were digitized and 
compared with each other. An average SPE pulse- 
shape was deduced, replacing the simple model that 
used only fixed rise- and decay-times. 

• The scintillation light was traced through the cell 
to the PMTs. SPE pulse-shapes with constant 
charge-to-amplitude ratio were added up to the fi- 
nal PMT pulse for each photo-electron produced in 
the cathode. The trigger threshold was compared 
to the amplitude of the total PMT pulse. However, 
the charge-to-amplitude ratio of the measured SPE 
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FIG. 2. Tlie upper plots sliow the simulated and measured 
trigger efficiency for low and high thresholds as a function of 
energy deposited in the center of one cell. Dots represent data 
while the solid line shows the simulated efficiency. The lower 
plots show the energy corresponding to a trigger efficiency of 
50% for each cell. The spread between data and Monte Carlo 
has been improved by a factor of about two compared to H. 



pulses varied slightly from pulse to pulse, result- 
ing in a smeared trigger threshold when plotted 
as a function of ADC counts. Instead of varying 
the width of the average SPE pulse in the simula- 
tion, the relative height of the threshold was sam- 
pled from a Gaussian distribution with fj, = 1 and 
a adjusted to describe the measured slope of the 
trigger-efficiency versus ADC-counts. The result- 
ing modeling improvements can be recognized by 
comparing the trigger threshold accuracy in Fig. g 
with the analogous quantities in Fig. 9 of 0). 

• The constant thresholds used for all cells and runs 
were replaced by an individual threshold for each 
discriminator. Variations in time were also taken 
into account by tracking thresholds using neutrino 
runs. 

In conclusion the new simulation reduced the spread of 
data/simulation over all cells in the detector (Fig. 0) from 
19.2% to 10.2% for the low (SPE) threshold and from 
7.6% to 3.5% for the high threshold. 



C. Test of the z/c detection efficiency 

^^Na and Am-Be sources were used to verify the abso- 
lute efficiency of tlie detector for positron annihilations 
and subsequent neutron captures. The 1.275 MeV pri- 
mary 7 of the ^^Na source is accompanied 90% of the time 
by a low energy positron which annihilates in the source 
capsule. The primary 7 mimics the positron ionization 
associated with a low energy v^ event and, together with 
the annihilation 7's, closely approximates the positron 
portion of a v^ event near the trigger threshold. 
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AmBe Efficiency at 43 Locations 
Data/MC: n=1.021 0=0.035 
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FIG. 3. Comparison of data (points) and Monte Carlo (his- 
tograms) for detection efficiency for ^^Na and Am-Be source 
runs at various locations. For positions of the radioactive 
source near the border of the central detector we measure 
lower efficiencies in good agreement with the simulation (see 
locations 3, 5, 22, ... for ^'^Na or locations II, 12, 18 ... for 
Am-Be). 

The ^^Na source was inserted into the central detec- 
tor at various locations during four dedicated calibration 
periods, separated by several months. A total of 85 differ- 
ent runs were taken in order to sample various distances 
from the PMTs and edges of the fiducial volume. This 
allows determination of an absolute efficiency since the 
source activity is known to 1.5%. After applying loose 
cuts to suppress background and correcting for detector 
deadtime, the measured absolute trigger efficiency could 
be compared with the Monte Carlo prediction; the results 
are shown in the top panel of Fig. |3|. Good agreement 
is seen in the average efficiency over all runs (the spread 
in data/Monte Carlo has been improved to 6.7% in 85 
locations compared to 11.1% in 36 locations in H), and 



the agreement between the four calibration periods was 
better than 1.4%. The ^^Na energy spectra predicted by 
the simulation and measured in the data also agree well. 
This comparison tests all aspects of the simulations: the 
high and low trigger thresholds, and the total energy 
deposit. 
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FIG. 4. Comparison of data (points) and Monte Carlo sim- 
ulation (histograms) for the spectra of total energy and first, 
second, and third most energetic hit (i?totai, Ei, E2, and Ea) 
for capture cosmic muon induced neutrons. The sharp fea- 
ture at 3.5 MeV total energy is related to the requirement 
that at least one sub-event has total energy above this value 
(see text). 

In order to check the neutron capture detection effi- 
ciency, the Am-Be neutron source was attached to one 
end of a thin (7.5mm) NaI(Tl)-detector readout by a flat 
PMT, so that the entire assembly could be still inserted 
in the gaps between each cell and its neighbor above or 
below. The NaI(Tl)-detector tagged the 4.4 MeV 7 emit- 
ted in coincidence with a neutron. 

The Nal(Tl) tag forced the digitization of the 4.4 MeV 
7 as the prompt part of an event and opened a 450 /iS 
window for neutron capture, the same coincidence win- 
dow as used in the i>e runs. 

Loose neutron cuts were applied, and corrections were 
made for detector deadtime and a low rate of random 
background. On average, the Monte Carlo efficiency pre- 
dictions agree well over the 43 locations tested (compared 
to 25 locations in Q ) with an average agreement of better 
than 2.1%, as shown in the lower panel of Fig. ^. Small 
adjustments of parameters of the detector simulations 
could improve the agreement of the AmBe efficiencies, 



but at the same time led to larger disagreement for other 
parameters, e.g. ^^Na efficiencies or the shape of energy 
spectra. Therefore, only directly measured parameters 
(trigger efficiency as a function of charge and shape of 
PMT pulses) were used to adjust Monte Carlo parame- 
ters. Measured efficiencies for ^^Na and AmBe were only 
used to estimate the systematic uncertainty of the simu- 
lation. 

Again, the energy spectra for neutrons predicted by 
the simulation and measured in the data were compared. 
The total energy seen in all cells and the energy detected 
in the three most energetic hits is plotted in Fig. ^ This 
test was done with cosmic muon induced neutrons, which 
are the dominant type of correlated events in neutrino 
runs. The neutrons were equally distributed over the 
whole detector. 

These procedures completely test our Ug efficiency sim- 
ulation. Thus, our ability to accurately generate the 
events, model the detector response, reconstruct the 
events, and correctly calculate the livetime of the data 
acquisition (DAQ) system was verified. 

The Monte Carlo simulation yielded an average effi- 
ciency over the whole detector as a function of i>c energy. 
The simulation included interactions in the acrylic walls 
of the cells, since there is significant efficiency for inverse 
beta decay originating there. Next, the efficiency from 
the simulation was folded with the incident v^ spectrum 
(including possible distortions due to oscillations) to get 
the overall efficiency which generally depends on the os- 
cillation parameters Am^ and sin^ 29. 



IV. ANALYSIS AND RESULTS 

We briefly discuss here event reconstruction, event se- 
lection, efficiencies, and backgrounds. Details may be 
found in [y||]. Both the analysis presented here and the 
one used for our previous papers were repeated without 
changing selection cuts for the present data-set. 

The energies and positions associated with hits were re- 
constructed for each bank. The position of the hit along 
the length of the cell was determined from TDC times 
with a time-walk correction applied on the basis of the 
collected charge. The collected charge for each end was 
corrected for light attenuation and PMT non-linearity 
and converted to an energy using energy calibration con- 
stants. The hit energy was determined as the weighted 
average of the measurements from either end. 

To select events in the energy ranges where the triggers 
are efficient, we required that each sub-event (prompt 
"positron" and delayed "neutron") have at least one hit 
greater than 1 MeV and at least two additional hits with 
energy greater than 30 keV. Any event with hits greater 
than 8 MeV in either sub-event was discarded. The mag- 
nitude and pattern of energy deposits in the prompt sub- 
event were required to resemble what was expected from 



the kinetic energy of the positron and its annihilation. 
(The annihilation 7's each had to have energy less than 
600 keV, and together less than 1.2 MeV. This is the only 
cut which treats the two sub-events asymmetrically.) The 
prompt and the delayed sub-events were required to be 
correlated in space and time. To further suppress back- 
grounds, an event was accepted if it started at least 150 
/iS after the last veto hit and at least 3.5 MeV of energy 
was deposited in either the prompt or delayed sub-event. 

The event yield must then be corrected for the effi- 
ciency of trigger and selection cuts as well as for detector 
deadtime, which has two components. The first one is 
the loss of neutrino events due to muons crossing the de- 
tector (a) within 150 jjls before the start of the neutrino 
event or (b) between the prompt and delayed sub-events. 
Its magnitude was determined from the measured muon 
veto rate and the distribution of inter-event times from 
detector simulation. The experiment livetime after losses 
due to the muon veto is approximately 66%. The second 
deadtime component comes from the DAQ system being 
unavailable to digitize a triple. Its magnitude is the ra- 
tio of the number of triples for which the DAQ was busy 
to the total number of triples seen by the trigger which 
could be directly measured using scalers. The deadtime 
of the trigger itself was measured to be less than 0.1%. 
The experimental DAQ livetime was about 81% for 1998 
and 92% for 1999-2000. The higher DAQ livetime in 
1999-2000 was due to recording correlated events only, 
rather than all triples, thus strongly reducing the load on 
the DAQ system. For the case of no oscillations, the com- 
bined efficiency of the trigger and selection cuts on neu- 
trino interactions is about 18%. The detector deadtime 
further reduces the efficiency to about 10% (the exact 
figure for each period being given in Table I). 

Experimental backgrounds may be naturally classified 
as uncorrelated and correlated^ with uncorrelated back- 
grounds due to random coincidences between triple trig- 
gers within the delayed coincidence window, and corre- 
lated background due to events in which both sub-events 
belong to the same process. 

The dominant source of uncorrelated events is natu- 
ral radioactivity. The inter-event time distribution for 
uncorrelated background events follows an exponential 
function with a time constant of ^500 /iS, as would be 
expected given the muon veto rate of ~2 kHz and the 
veto-dependent event selection requirements. This time 
dependence is slow compared to that of signal and cor- 
related backgrounds, hence the uncorrelated background 
could be separated and studied by looking at long inter- 
event times. 

The main source of correlated background are neu- 
trons from muon spallation or capture. These events are 
mainly comprised of proton-neutron events in which a 
single neutron deposits its kinetic energy by scattering 
from protons and is then captured, and double neutron 
events in which two (typically thermal) neutrons from 



TABLE I. Data taking periods, efficiencies (including livetime), measured event rates A''i and N2, and results of the swap 
analysis (see text), including the various background estimates. Uncertainties are statistical only. 



Period 


1998 


1999-1 


1999-11 


2000 


Reactor 


on 


890 m off 


on 


750 m off 


on 


890 m off 


on 


890 m off 


time (days) 


30.4 


29.4 


68.2 


21.8 


60.4 


29.6 


83.2 


27.5 


efficiency (%) 


8.0 


8.0 


11.5 


11.6 


11.6 


11.6 


10.9 


10.8 


measured rates 


iVi (d-^) 


39.6 ± 1.1 


34.8 ± 1.1 


54.9 ± 0.9 


45.1 ± 1.4 


54.2 ± 0.9 


49.4 ± 1.3 


52.9 ± 0.8 


43.1 ± 1.3 


N2 (d-i) 


25.1 ± 0.9 


21.8 ± 0.9 


33.4 ± 0.7 


32.0 ± 1.2 


32.5 ± 0.7 


32.6 ± 1.0 


30.2 ± 0.6 


30.4 ± 1.1 


(l-ei)Bpn (d-i) 


0.88 


0.89 


1.11 


1.11 


1.11 


1.11 


1.07 


1.07 


efficiency corrected rates 


Background (d"'^) 


292 ± 11 


255 ± 10 


265 ± 6 


266 ± 10 


256 ± 6 


265 ± 9 


249 ± 5 


272 ± 9 


R. (d-i) 


202 ± 19 


182 ± 18 


212 ± 10 


124 ± 17 


214 ± 11 


161 ± 15 


237 ± 10 


129 ± 16 


i?calc (d-^) 


216 


154 


218 


129 


220 


155 


218 


154 



the same spallation event are captured in the detector. 

Our analysis is based on 350.5 days of data taking, of 
which 242.2 days were at full power and the remainder 
at partial power with a reactor down for refueling. For 
the analysis we subdivided the data into 8 periods. Four 
of the periods correspond to the four reactor refueling 
periods in which one reactor was off (o/f periods). Each of 
the remaining four periods {on periods) are constructed 
from intervals of full power data bracketing each refueling 
period. Table B shows the running time for each of the 8 
periods and the distance to the down reactor for each of 
the Oj(f periods. 

The raw trigger rates for triples and correlated triples 
were approximately 50 Hz and 1 Hz, respectively. For 
1999-2000, the typical event rate after selection was '^55 
d""'^ with all reactors at full power. Under the assump- 
tion of no oscillations, the efficiency after the trigger, 
deadtime, and event selection for detecting v^s above in- 
verse beta decay threshold was ~11%; precise estimates 
of the efficiency period-by-period are listed in Table |. 
The observed event rate, TVi, may be compared to an 
expected signal rate of '^25 d^^ for no oscillations, im- 
plying a signal-to-noise ratio of ^ 0.8. The uncorrelated 
background event rate after selection was ^1 d~^. 



A. Analysis Avith the "reactor power" method 

From Table it is evident that the event rate is signif- 
icantly lower during each refueling period. 

To investigate more quantitatively the correlation be- 
tween event rates and reactor power, we plot in Fig. ||the 
experimental rate corrected for efficiency and deadtime 
-Roxp, against the calculated signal rate i?caic expected 
for no oscillations. Only statistical uncertainties are indi- 
cated. If the data were consistent with no oscillations and 
the background were constant, then the points should lie 
along a straight line with unity slope. The y-intercept 
is equal to the rate of background interactions scaled by 
the ratio of the effective background detection efficiency 



to the neutrino detection efficiency. The data are in fact 
consistent with lying along a straight line. A linear fit to 
these data gives a slope of 1.011 ± 0.104 (stat.) and a 
y-intercept of 257.5 ± 20.7 (stat.) d"^ The reduced x^ 
of the fit is 0.89. Our data are therefore consistent with 
the hypothesis of no oscillations. 
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FIG. 5. The event rates i?cxp for different data taking peri- 
ods, corrected for deadtime and neutrino detection efficiency, 
plotted versus the expected neutrino interaction rate 7?caic 
for no oscillations. Errors are statistical only. Points cor- 
responding to data taking periods with same reactor power 
conditions should lie on top of each other. Also shown is the 
result, discussed in the text, of a linear fit to the data. 

We have also analyzed the energy dependence of the 
neutrino interactions to see whether it is consistent with 
no oscillations as well. For each of the four pairs of on- 
off periods, one may subtract the event rate at partial 
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FIG. 6. The prompt energy spectrum after on- Oj(f subtrac- 
tion averaged over the 4 pairs of on/ off periods. The his- 
tograms show the corresponding expectations for no oscilla- 
tions (solid line) and the Kamiokande best fit (dashed line). 



power from the rate at full power. The resulting differ- 
ence, after the small correction for the fuel burn-up has 
been made, is the contribution to the full power event 
rate from the neutrinos emitted by the reactor unit which 
was being refueled during the off period. Fig. y shows 
the measured on-off event rate difference binned in the 
visible prompt (positron) energy. These data are not cor- 
rected for efficiency or deadtime. To make this plot, the 
weighted average of the four pairs of on-off periods is 
taken. Also shown are the corresponding expectations 
from Monte Carlo for two scenarios: no oscillations and 
oscillations with parameters obtained from the best fit 
to the Kamiokande data |^] . The comparison of our data 
with Monte Carlo shows that the hypothesis of no os- 
cillations is consistent not only with the measured event 
rate but also with the observed positron energy spectrum 
(x^/d.o.f = 1.39 for 8 degrees of freedom), while it is not 
consistent with the rates implied by the Kamiokande best 
fit parameters (x^/d.o.f = 3.69). 

To test our data for oscillation hypotheses throughout 
the Am^ — sin 29 plane for two flavor mixing, a x^ ana- 
lysis using the "reactor power" changes was carried out. 
The x^ is defined as 
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where -Rexp is the observed rate for period j, bg is the 



background rate, i?caic is the calculated rate for the pe- 
riod i which depends on Arri^ and sin 26, and a accounts 
for possible global normalization effects due to systematic 
uncertainties; Usyst is 0.061 (see Table O and the section 
on systematic errors below) . The quantity bg is scaled by 
0.94 for the two periods in 1998 to account for the dif- 
ferent trigger conditions in 1998 but is otherwise treated 
as constant. This scaling factor was determined from 
measuring how the rates for double neutron background 
events (the selection of which is described in the follow- 
ing section) and uncorrelated background events changed 
between 1998 and 1999. 

To define the 90% confidence level (CL) acceptance 
region for our data, we have followed the procedure sug- 
gested by Feldman and Cousins (FC) |g^. We imple- 
mented this procedure in two ways. 

The first, called Monte Carlo implementation, was re- 
alized as follows. A Xbcst ^^s determined by minimiz- 
ing Eq. (||) with respect to bg, a, Am^, and sin^ 20 
for physically allowed values of Atti^ and sin^ 29. We 
found the best fit to correspond to a null sin^ 29 with 
a slightly un-physical a = 1.007 and a xtcst/^-'^-f- ~ 
5.8/6 . The physical Am^ — sin^ 29 plane was then sub- 
divided into a fine grid. At each grid point, we min- 
imized Eq. (n[) with respect to bg and a to obtain 
AXdata = Xlata(^"^^sin2 26') - xLsf For determin- 
ing whether the grid point was allowed at the 90% CL, 
we simulated IC* independent experiments at each grid 
point. The same x^ minimization procedure was car- 
ried out for each simulated experiment as for the data 
to obtain 10'* Ax^c's- These Axle's were sorted in in- 
creasing order to find Axc , the value of Ax^ greater than 
90% of the Axmc's- ^^ ^Xdata < ^xli the grid point was 
accepted. 

The second way in which we implemented the FC 
procedure, called raster scan, subdivides the two- 
dimensional grid in Arri^ slices. For each value of 
Am'^, Eq. (|l|) is minimized with respect to bg, a, and 
sin^ 26. The value obtained for sin^ 29, without restrict- 
ing the fit to the physically allowed range, is denoted 
(sin" 20)bost a-iid its corresponding error is denoted asm- 
The best fit is slightly unphysical; for Arn^ = 0.1 eV^, 
(sin" 26)hcst I o' sin = —0.2. The one-dimensional 90% CL 
upper limit on sin 29 at the fixed value of Ani^ is then 
given by apco'im where apc is looked up from Table X 
of Reference |27) for xq = "^ °* . 

While the raster scan method does not yield the global 
X^ minimum in the Am^ — sin 29 plane, it is computa- 
tionally much faster. Checks have been carried out that 
the two methods for implementing the Feldman-Cousins 
procedure yield the same limits. For the purpose of de- 
termining the regions of parameter space excluded by 
our data, knowledge of the x^ global minimum is not re- 
quired. We have therefore used the raster scan method 
to obtain the exclusion curves reported in this paper. 



The dashed curve in Fig. shows the region of /S.iv? — 
s\T? 26 plane excluded at the 90% CL by our data an- 
alyzed with the "reactor power" method. In the limit 
of large Am^, the range sin^ 26* > 0.33 is excluded; 
whereas in the limit of maximal mixing, the range Arri^ > 
1.6 X 10^'^ eV^ is excluded. We note that, in the limit of 
large Arn^, the Monte Carlo method excludes the range 
sin^ 261 > 0.35. 

As already mentioned the independent analysis dis- 
cussed in detail in p^pl was also improved and repeated 
for the full data-set. In this case the "reactor power" 
analysis differs from the one described above in that the 
data are more finely binned by run rather than averaged 
by period. There were typically two runs per day. Af- 
ter combining short runs (runs with fewer than 6 neu- 
trino candidates) with adjacent runs, 698 data points 
were obtained. A x^ analysis identical in approach to 
that described above was carried out using a system- 
atic error that in this case amounts to 6.9%. Again, 
the best fit is slightly unphysical; for Am^ = 0.1 eV^ 
(sin^20)bcst/o's,;„ — —0.5. The 90% CL exclusion con- 
tour obtained in this analysis is very similar to the dashed 
curve in Fig. 0, but it is shifted toward smaller sin^ 29, 
with sin^ 26 > 0.29 excluded in the large Am^ limit. The 
shift in the exclusion boundary is consistent with small 
systematic differences expected between the two indepen- 
dent reconstructions and analyses. 



B. Analysis with the "swap" method 

The "swap" method, where the background is directly 
subtracted rather than using modulation of the reactor 
power, has substantially greater statistical power than 
the "reactor power" and, in addition, it has somewhat 
different systematics. 

We briefly describe the method here; detailed descrip- 
tions have already been published [p],p|j2q|. Let A^i be 
the event rate after applying the neutrino selection cuts 
described above. We then call N2 the rate obtained by 
applying the positron cuts to the delayed sub-events and 
the neutron cuts to the prompt sub-events ("swapped" 
selection) . The measurements of A^i and N2 are listed for 
each period in Table |. It is found that only about 20% 
of the neutrino signal cancels in the difference Ni — N2, 
as determined from Monte Carlo simulation. At the 
same time the uncorrelated background and the double 
neutron component that dominates the correlated back- 
ground cancel in the difference. We call (1 — ei)i?pn = 
(1 - esp)Bpn,sp + (1 - ecap)Spnxap the rcsidual contri- 
bution to A^i — A^2, mainly due to the proton-neutron 
("pn") component of the correlated background. Here 
the e's refer to the efficiency of the "swapped" selection 
for each channel pq], "sp" denotes neutron production 
by fi spallation (mainly in the laboratory walls), "cap" 



by capture (mainly in the water buffer) and "1", main- 
taining the notation from our earlier papers, the total. 
While the latter process is well understood and can be 
reliably calculated with Monte Carlo simulations, the for- 
mer is rather poorly known. The shape of the prompt en- 
ergy spectrum for neutrons from spallation was obtained 
generating neutrons in the laboratory walls according to 
several parametric models and passing them through the 
detector simulation and event selection. The set of para- 
metric models spanned the range of uncertainty in our 
knowledge of the energy dependence of neutron produc- 
tion. The normalization was then determined by assum- 
ing that high-energy neutrino-like events, selected by re- 
placing the cut on maximum hit energy (< 8 MeV) in the 
"positron" sub-event by the requirement that at least one 
prompt hit have energy greater than 10 MeV, are due 
exclusively to spallation. The uncertainty on the energy 
spectrum, quantified by the dispersion between the dif- 
ferent models, was taken into account in the systematic 
error. The term -Bpn.cap was found from the measured 
muon rate through the detector combined with the veto 
inefficiencies and the relatively well known total neutron 
production cross section and energy distribution by muon 
capture on oxygen. The main systematic derives from the 
veto inefficiency. 

The quantity (1 — ei)i?pn was estimated period by pe- 
riod, and the results are shown in Table ffl. The mag- 
nitude of (1 — ei)i3pn is small compared to the differ- 
ence Ni — N2, that is, the contribution from the proton- 
neutron component of the background largely cancels in 
the difference. Therefore, even though the systematic 
error on (1 — ei)i?pn is of order 100%, the resulting con- 
tribution to the systematic error on the neutrino signal 
is only a few percent. In the sixth line of Table ffl we list 
for each run-period the resulting background (assuming 
for the purpose of this illustration that the background 
efficiency is the same as for the signal, and correcting for 
livetime). The observed Dc rate (i?i/), corrected for live- 
time and efficiency, and the expected neutrino rate i?caic 
for no oscillations are also given in the Table. 

Similarly to the reactor power analysis, we have carried 
out a x^ analysis to test our data for oscillation hypothe- 
ses throughout the two flavor oscillation Am^ — sin^ 26 
plane. The x^ definition is 
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where Usyst for the "swap" method is estimated to be 
0.053 as discussed below. The free parameters in this 
definition of the x^ are Am^, sin^ 26, and a. The Monte 
Carlo method gives xlest/'^-^-f- = 10-3/7 for sin^ 20 
consistent with zero and, again, a slightly un-physical 



a = 1.008. 

The region of parameter space excluded at the 90% 
CL by this analysis, based on the raster scan method, is 
indicated by the sohd curve in Fig. |7[ In the hmit of large 
Ato^, the range sin^ 26* > 0.164 is excluded; whereas in 
the limit of large mixing, the range Am^ > 1.1 x 10^"^ cV'^ 
is excluded. We note that, in the limit of large Ato^, the 
Monte Carlo method excludes the range sin^ 20 > 0.162, 
and gives an essentially identical exclusion curve. 
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FIG. 7. Regions of Am^ — sin^ 2Q plane (two flavor oscilla- 
tions) excluded at the 90% CL by the "reactor power" analysis 
(dashed curve) and "swap" analysis (solid curve). Also shown 
are the Kamiokande allowed region and best fit (star) and the 
region excluded by the Chooz experiment 



C. Test of the "swap" method 

As a further test of the "swap" method, we have inves- 
tigated the energy dependence of A^i and N^- The mea- 
sured energy dependence was compared to what would 
be expected on the basis of our assumptions about the 
signal and background. Were a significant source of back- 
ground ignored or incorrectly treated, a discrepancy be- 
tween data and expectation would result. 

To carry out this investigation, we assembled five sam- 
ples of events: 

• Vc- Inverse beta decay events were generated and 
simulated in the detector with normalization deter- 
mined from the reactor powers, cross section, and 



number of target protons. No oscillations were as- 
sumed, as consistent with the outcome of the reac- 
tor power analysis. 

• uncorrelated background: These events were se- 
lected from our data by inverting the spatial 
and temporal correlation requirements between the 
prompt and delayed sub-events. The data sample 
was normalized to reproduce the event rate at large 
inter-event times. 

• B-p-cis-p'- Neutrons produced by muon spallation in 
the laboratory walls were generated and passed 
through our detector simulation. As already men- 
tioned, the sample was normalized by assuming 
that high-energy events satisfying the neutrino se- 
lection cuts are due to spallation. 

• .Bpn^cap^ Neutrons produced by muon capture in 
water were generated and passed through the detec- 
tor simulation. This data sample was normalized 
on the basis of the measured muon rate through the 
detector, the muon veto inefficiency, the fraction of 
muons stopping in the water, and the cross section 
for muon capture and neutron emission. 

• double neutron: Double neutron events were se- 
lected from data by requiring a muon hit within 
100 ^s preceding the start of the event and applying 
the neutron capture cuts to both the prompt and 
delayed sub-events. The sample was normalized 
so that-after application of the neutrino selection 
cuts-the combined 5 samples gave the measured 
total A^i rate. 

The five data samples were subjected to the neutrino 
event selection cuts (A^i) and the swap event selection 
cuts (A'^2), respectively, and summed. The resulting en- 
ergy spectra, with statistical uncertainties, are shown as 
histograms in Fig. y. The expectation from the sum 
of the five samples is in good agreement with the data 
(points). Keeping in mind that only the overall normal- 
ization of the A^i spectrum is not independently deter- 
mined, the validity of the "swap" analysis is solidly sup- 
ported by this test. 



D. Systematic uncertainties 

The systematic uncertainty receives contributions from 
the detection efficiency and the flux calculation. In ad- 
dition, the "reactor power method" suffers a systematic 
error from background variations with time, and the sys- 
tematic uncertainty in the "swap" method has a contri- 
bution from the uncertainty in the estimate of (1— ei)i?pn. 

We have estimated the systematic uncertainty in the 
detection efficiency as follows: 







, ^ Nj (normal cuts) 






• 1 — 1 
• 


• Data 


- 




• 


« — Expected: V + backgr. 




"^ 


• 



8 9 10 

Energy (MeV) 




9 10 

Energy (MeV) 

FIG. 8. The energy spectrum of events comprising TVi and 
A^2- The points are measurements while the histogram shows 
the result of calculations described in the text. The shading 
indicates the statistical uncertainty in the calculation, while 
error bars for the measurements are too small to be visible. 



• i>o selection cuts efficiency: The neutrino event se- 
lection cuts were varied randomly in the multi- 
dimensional cut space over a reasonable range. For 
each variation, the ratio of the observed number of 
events to that expected for no oscillations was cal- 
culated. The systematic error in the event selection 
efficiency was taken to be the rms of the variations 
in the ratio. The uncertainties in the definition of 
the energy scale are absorbed in this error com- 
ponent. The lower systematic uncertainty for the 
swap method is due to cancellation of some system- 
atics in the difference iVi ~ N2. 

• e"*" trigger efficiency: The systematic uncertainty 
in the e"'" trigger efficiency is based on comparison 
of simulated efficiencies with the measured efficien- 
cies for the ^^Na calibration runs described above. 
To decouple uncertainties in the event selection ef- 
ficiency from uncertainties in the trigger efhciency, 
loose cuts designed to have negligible inefficiency 
were applied to select ^^Na events for this analy- 
sis. The run-by-run comparison of the simulated 
and measured efficiencies has already been shown 
in Fig. H Averaged over all runs, the efhcicncics 
agree to 0.2% with a rms of 6.7%. Grouping the 
runs by four calibration periods, the agreement was 
1.3%. Combining this with the 1.5% uncertainty in 
the activity of the source, the estimated systematic 



error is 2.0%. 

• n trigger efficiency: In an approach similar to that 
for estimating the systematic error in the e"*" trigger 
efficiency, we have used the Am-Be calibration runs 
described above to estimate the systematic uncer- 
tainty in the n trigger efRciency The run-by-run 
comparison of measured versus simulated n trig- 
ger efficiency has already been shown in Fig. y. 
The simulated efficiency is typically lower than the 
measured one. This difference is largely systematic, 
as manifested by the relatively small rms of 3.5% 
across the different calibration periods and posi- 
tions. Averaged over runs, the difference between 
simulation and measurement is 2.1%, which we as- 
sign as the systematic uncertainty for the n trigger 
efficiency. 

The results for these sources of systematic error for the 
two analysis methods are listed in Table |l[ As explained 
in Section II. B the systematic uncertainty in the Dc flux 
is estimated to be 2.1%. 

The stability of background rates is a key assumption 
for the "reactor power" analysis. The actual level of back- 
ground stability was estimated by comparing the average 
rate (livetime-corrected) during the full power periods 
to the average rate during the partial power periods for 
several background data samples: double neutrons, Bp^, 
and uncorrelated background. 

In addition, Michel electron events, present in the data 
due to the inefficiency of the veto detector, were used to 
track changes in the veto efficiency and, in particular, 
in the background due to neutron production by muon 
capture in water. These events are selected by requiring 
no activity in the muon detector, energy depositions in 
the prompt sub-event consistent with a muon track, and 
a delayed sub-event 5-20 fis later with an energy deposit 
of 10-70 MeV. 

The Michel electron data sample was observed to have 
a rate stability better than 5% and all other data sam- 
ples were found to be stable to better than 1%. The 
rate variation for each background was normalized to 
its estimated contribution to the neutrino event rate 
and then divided by 12.9 d~^, the average difference in 
livetime-corrected event rates between full power and 
partial power periods. Combined, the resulting ratios for 
the four backgrounds indicated a background instability 
of 2.1% relative to the signal. We this take 2.1% as our 
estimate of the contribution to the systematic error in the 
"reactor power" analysis from background variations. 

The contribution to the uncertainty in (1 — esp)^pn,sp 
from muon spallation in the walls was estimated from 
the spread in results from the four different models used 
to simulate neutron production and was found to be 
0.29 d^^ . The estimated contribution to the uncertainty 
in (1 — ecap)^pn,cap from the veto counter inefficiency re- 
sulted to be 0.94 d~^. Thus, the total systematic uncer- 
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TABLE II. Contributions to the systematic error of the 
"reactor power" and "swap" analyses. 



Error source 


"reactor power" 


(%) 


"swap" (%) 


e^ trigger efficiency 


2.0 




2.0 


n trigger efficiency 


2.1 




2.1 


Uc flux prediction 


2.1 




2.1 


Uc selection cuts 


4.5 




2.1 


Background variation 


2.1 




N/A 


(1 — ei)-Bpn estimate 


N/A 




3.3 


Total 


6.1 




5.3 
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tainty on (1 — ei)i?pn amounted to 0.98 d~^. This result 
was livetime-corrected and corresponded to 3.3% of the 
average livetime-corrected value for A^i — N2- 

The individual contributions are shown in Table O and 
added in quadrature to obtain the total systematic error 
for each analysis method. 



V. CONCLUSION 

The results presented here, based on nearly double the 
number of events of our previously published Palo Verde 
data, confirm the absence of v^ — > v^ oscillations for low 
energy reactor neutrinos. The excluded regions for our 
"reactor power" and "swap" methods are enlarged ac- 
cordingly. For the "reactor power" method the new mix- 
ing angle limit is only slightly more restrictive than in 
our previous results. This is due to a small shift of the 
central value of the fit. A substantially larger exclusion 
region is obtained with the "swap" method, thanks to 
reduced systematics. In conclusion we find that the ratio 
of observed interaction rate to the one expected for no os- 
cillations is i?obs/i?caic = 1.01±0.024(stat)±0.053(syst). 
These final results are dominated by systematics errors. 

Our measurements, along with those reported by 
Chooz P,p| and Super-Kamiokande H, excludes two 
family v^ — v^ mixing as being responsible for the atmo- 
spheric neutrino anomaly reported by Kamiokande B. 
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